Slabs of stabilized jellium: Quantum-size and self-compression effects 
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We examine thin films of two simple metals (aluminum and lithium) in the stabilized jellium model, 
a modification of the regular jellium model in which a constant potential is added inside the metal 
to stabilize the system for a given background density. We investigate quantum-size effects on the 
surface energy and the work function. For a given film thickness we also evaluate the density yielding 
energy stability, which is found to be slightly higher than the equilibrium density of the bulk system 
and to approach this value in the limit of thick slabs. A comparison of our self-consistent calculations 
with the predictions of the liquid-drop model shows the validity of this model. 
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I. INTRODUCTION 

Thin films or slabs are systems made out of a few lay- 
ers of atoms: they are finite in one direction and infinite 
in the other two perpendicular directions. Various the- 
oretical models are available to calculate the electronic 
structure of slabs. One of the important features coming 
out of these-|Calculations is the so-called quantum-size ef- 
fect (QSE),El i.e., the influence of the finite size on various 
physical properties of the slak These effects, which can 
be experimentally recognizedjj decrease as the size of the 
slab increases. In fact, surface energies and work func- 
tions of the semi-infinite system are often derived from 
thin-slab calculations, which are simply extrapolated to 
this limit .ETeI 

The simplest model to predict the electronic structure 
of simple sp-bonded metals is the jellium model, where 
the ions are replaced by a positive neutralizing back- 
ground. Within thiSppodel, the QSE of thin films was 
examined by SchulteE He found an oscillatory behavior 
of the work function, as a function of the thickness of 
the slab. The same oscillatory behavior is found for the 
surface energy, defined as the energy required, per unit 
area of aew surface formed, to split the solid in two along 
a plane .0 

The jellium model has been referred to as giving-, ip^ 
sight into the realistic QSE appearing in real systems.ETtS 
Notwithstanding important differences, an oscillatory 
pattern also appears in atomistic first-principles slab cal- 
culations of both the work function and the surface en- 
ergy. Although the presence of the lattice may obscure 
the periodicity and the amplitude of the QSE, extrema 
were found at positions which agree with the jellium re- 
sults. However, in the case of first-principles calculations 
difficulties arise (due to the cumbersome numerics) when 
one is to extract well- converged surface properties from 



thin films made of typically 2 to 15 layers.ErQ Hence, the 
clean jellium QSE, with no uncertainties in the extrapo- 
lated results, remains as a guide for more realistic inves- 
tigations. 

In this paper, we consider slabs in the framework of 
a simple modification of the jellium model which yields 
energy stability against changes in the_background den- 
sity. This so-called stabiHzed jelhumEil or structureless 
pseudopotential model yields reahstic results, especially 
in the case of metals with high valence-electron density. 
For instance, the stabilized jellium model predicts pos- 
itive surface energies that increase rapidly at high elec- 
tron densities, as shown by experiment, while the jellium 
model predicts surface energies that are strongly neg- 
ative at these densities. The stabilized jellAwn model, 
first introduced by Perdew, Tran, and Smithll3 and sim- 
ilar to-|tlie ideal-metal concept developed by Shorej-SLad. 
Rose,t3'til hajS-jbeen applied to the study of surfaceala'tSi 
and clusters.EZI In a way, the stabilized jellium is in be- 
tween the jellium model and more sophisticated atomistic 
approaches: although it is still a continuous model (one 
may choose slabs with arbitrary thickness) with an an- 
alytical expression for the bulk energy, its physical pre- 
dictions are in reasonable agreement with experiment. 
Besides including electrostatic corrections to the jellium 
model, the stabilized jellium model contains an averaged 
pseudopotential correction. 

We calculate the self-consistent energetics (surface en- 
ergy and work function) of slabs of stabilized jellium, 
with use of the local-density appfpximation (LDA) of 
density- functional theory (DFT).t3'llj We take two met- 
als: Al {r, = 2.07, Z 3) and Li {r, = 3.24, Z = 1), 
investigate the QSE, and compare our self-consistent slab 
calculations with those obtained for a semi-infinke stabi- 
lized jellium. We also test an extrapolation rule,c2l which 
has already been used to describe non-local surface ener- 
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gies of the bounded electron gas.0 

Although the stabilized jellium niadei j-ean be tai- 
lored to give face- dependent results, ll3il3'E2l it cannot 
describe the inhomogeneous relaxation predicted by 
first-principles calculations where the distances between 
atomic planes of the same family are optimized. How- 
ever, an interesting effect displayed by the stabilized jel- 
lium model, which cannot be accounted for Jay the jel- 
lium model, is the so-called self-compressinnrj (or self- 
expansion, in the case of charged systemscJo) of clus- 
ters. This effect, which can be classically viewed as the 
compression of a finite system due to the surface tension, 
is most prominent for small systems and almost negligible 
in the case of large clusters. 

We investigate here the self-compression of thin films. 
We fix the size of the system along the direction perpen- 
dicular to the surface, and search for the background den- 
sity which minimizes the total energy per valence electron 
of the slab. The equilibrium density is found to increase 
as the thickness of the slab decreases, and to converge to 
the bulk electron density in the infinite-thickness limit. 
Furthermore, the equilibrium electron-density parameter 
r* is found to oscillate with the slab thickness, as a mani- 
festation of the QSE, but the general trend is found t»]lp€j 
well described within the liquid-drop model (LDM)Ea'EZl 
based only on the knowledge of the bulk energy per unit 
volume and the surface energy. We discuss the relation- 
ship between this self-compression effect and the relax- 
ation of metal slabs predicted by atomistic first-principles 
calculations. 

In section II we present briefly the stabilized jellium 
model for slabs. In section III we discuss the results we 
obtained within this model. The main conclusions are 
drawn in section IV, where further comments on the re- 
lationship between the stabilized jellium and more elabo- 
rated models are made. Equations are written in atomic 
units throughout, i.e., = 7i = me = I. 



II. SLABS OF STABILIZED JELLIUM 

The stabilized jellium modclEl takes into account the 
lattice ions, but keeps the essential simplicity of the jel- 
lium model. The total energy is obtained as a functional 
of the electron density n(r), in the following way: 

Esj[n, 72+] = Ej[n, n+] + (cm + wr) J (fr n+(r) 

+ {Sv)^s J ^ ["W - n+{v)] , (1) 

where 

n+=ne(r) (2) 

represents a positive neutralizing background density, 
0(r) being a function which equals 1 inside a given sur- 
face and outside, and 
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(3) 



is the average valence-electron density. Ej is the regular- 
jellium total energy, cm is the Madelung energy arising 
from the Coulomb interaction between a uniform neg- 
ative background inside the spherical Wigner-Seitz cell 
and a point ion at its center. 



CM = 



10 r, ' 



(4) 



wn is the average value of the repulsive non-Coulomb 
part of the Ashcroft empty-core pseudopotential. 



(5) 



and {Sv)ws represents the difference between the local 
pseudopotential and the jellium potential, averaged over 
the Wigner-Seitz cell. 
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The core radius Vf. of the Ashcroft empty-core pseudopo- 
tential is chosen to stabilize the metal for given values of 
the electron-density parameter and the chemical va- 
lence Z. 

The two terms added to the regular-jellium energy Ej 
are a volume term and a surface term. They simply ac- 
count for the subtraction of the spurious self-interaction 
of the positive jellium background and the inclusion of a 
constant structureless potential inside the metal. This 
procedure may be understood as a first-order pertur- 
bation to a jellium system, but with the perturbation 
treated in an averaged manner. 

The density functional of Eq. (|]) represents the total 
energy of an arbitrary inhomogeneous system. In the case 
of an infinite uniform system, the equilibrium density is 
obtained from the bulk stability condition 



dr. 



= 0, 



where 



,bulk 



J}ulk 



ej + CM + WR 



(7) 



(8) 



represents the average bulk energy per valence electron, 
ej being the regular-jellium contribution. Within this 
model any individual metal minimizes the energy at 
a given equilibrium density, while the jellium energy 
presents a single minimum at 4.2 close to the 

electron-density parameter of sodium. 

We consider slabs of stabilized jellium. Slabs are trans- 
lationally invariant in the plane of the surface, which 
is assumed to be perpendicular to the z axis. Hence, 
the single-particle wave functions can be separated into 
a plane wave along the surface and a component ^(z) 
describing motion normal to the surface with energy e. 
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This component is obtained by solving self-consistently 
the Kohn-Sham equation 



'2dz 



2 , VHiz) + V^a{z) + Vp,{z) 



(z) = e0(z), (9) 



where Vh{z) represents the Hartree electrostatic poten- 
tial, Vxc{z) is the exchange-correlation potential, and 
Vps{z) accounts for the pseudopotential, 



(10) 



Vxc{z) is obtained in the LDA, using the^electron-gas cor- 
relation energy of Ceperl«2-.and Alder,E3 as parameter- 
ized by Perdew and Wang.c^ Essentially the same results 
are obtained-from the parameterizations p£ Vosko, Wilk, 
and Nusaiiu and of Perdew and Zunger.tll We have not 
chosen to use extensions-such as the generalized gradient 
approximation (GGA),l3 since the LDA has been shown 
to give surprisingly good results_in describing the prop- 
erties of jellium planar surfaces O 

Outside the positive background the electron-density 
profile n{z) decays rapidly from its bulk value n. The 
electronic system can therefore be taken to be finite in 
the z direction by assuming that n{z) actually vanishes at 
a given distance zg from the surface. Hence, we introduce 
infinite potential walls at a distance zq from each surface, 
and follow Ref. |^ to expand the wave functions (f>{z) in 
a Fourier sine series. The distance zq (typically 2 or 3 
Fermi wavelengths) and the number of sines kept in the 
expansion of the wave functions (t>{z) have been chosen to 
be sufhciently large for our calculations to be insensitive 
to the precise values employed. These calculations have 
been compared with others which we have carried out 
for a semi-infinite electron system by using the Monnier- 
Perdew codeE3 for the numerical integration of Eq. (^ . 

For a given thickness L of the slab, we obtain the sur- 
face energy from the difference between the total energy 
of Eq. (^) and the corresponding result for a homoge- 
neous electron gas of density n^, i.e.. 
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(11) 



where A is the normalization area. The work function is 
obtained as the difference between the computed values 
for the vacuum and Fermi levels of our electron system. 



III. RESULTS AND DISCUSSION 

First of all, we compare jellium and stabilized-jellium 
electron densities n{z) and effective potentials, 



Veffiz) = Vniz) + Vxciz) + Vpsiz). 



(12) 



Jellium and stabilized-jellium valence-electron densities 
and effective potentials for an Al slab of L = 2 Xp 

1 /3 

[Xp = (32 7r^/9) rs is the Fermi wavelength] are shown 



in Fig. 1, together with the positive background density 
n+. We note that the stabilized-jellium electron density 
is steeper at the two surfaces, so that the electronic spill- 
out is slightly smaller within this model. This is due to 
the fact that electrons feel a deeper effective potential. 
Both jellium and stabilized-jellium electron densities ex- 
hibit quantum oscillations inside the metal, the so-called 
Friedel oscillations J3 and an exponential decay outside. 

Figs. 2 and 3 show our calculated stabilized-jellium 
surface energies for slabs of Al and Li, respectively, as 
obtained from Eq. (11) versus the thickness L of the 
slab. Both curves show damped oscillations with minima 
occurring at the slab width L ^ nXp/^ (n = 1,2,...). 
The same QSE, which reflects the quantization of the 
electronic motion along cme direction, is known to occur 
within the jellium modeLu 

Both the average bulk energy per valence electron e'gj'^ 
and the surface energy of the semi-infinite stabilized jel- 
lium 



a = lim cr(L) 



(13) 



may be obtained from a linear fit of the following equa- 
tion: 



EsjjL) 
A 



= 2a + nLe' 



bulk 

SJ •> 



(14) 



where Es,j{L) represents the total energy of Eq. (|^). 
Following this procedure, we reproduce the bulk energy 
of Eq. (^) and predict surface energies of 925 erg/cm^ and 
311 erg/cm'^ for Al and Li, respectively. These surface 
energies, represented in Figs. 2 and 3 by horizontal solid 
lines, agree with those reported in Ref. ^for semi- infinite 
media. 

An alternative procedure to extrapolate the surface en- 
ergy a of the semi-infinite medium from our calculated 
thin- film surface energies (j{L) is to use the relatione^ 

»■= ^ , (15) 

where L„ represents the threshold width for which the 
nth subband for the z motion is first occupied. Analyt- 
ical insight for this procedure is encountered within the 
infinite-barrier model (IBM), where the effective poten- 
tial Veff{z) is replaced by an infinite square well and the 
one-particle wave functions 4>{z) are simply sines. Based 
on this procedure, the numerical error introduced in a 
by our slab calculations is found to be within 0.1%. The 
advantage of this algorithm is that we simply need three 
points to obtain the asymptotic limit, while the linear 
fitting may yield erroneous results if one only takes a few 
thin films. 

Slabs with L < 0.5 Aj^ are interesting in their own, 
since they can be constructed in the laboratory, e.g., by 
joining two different semiconductors. Nevertheless, we 
do not give results for these ultra-thin slabs, since they 
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fall within the two-dimensional limit where the three- 
dimensional LDA and GGA .farmulae for exchange and 
correlation are known to fail.E^ 

For comparison, first-principles thin-film calculations 
of the surface energy of the most dense faces of Al and 
Li [(111) for fee Al and (0001) for hep Li] are represented 
in Figs. 2 and 3 by solid circles and triangles, with the 
slab width of a z^-layer unrelaxed crystalline film taken 
to be V times the interplanar distance. For Al there is 
reasonable agreement between our stabilized-jellium re- 
sults and atomistic first-principles calculations, the am- 
plitude of the stabilized-jellium oscillations being com- 
parable to that exhibited by first-principles calculations. 
For Li, however, there is a serious discrepancy between 
stabilized-jellium and first-principles calculations. Since 
lithium has been found to behave to some extent like 
a covalent solid rather than a free-electron gasJ^ZTcy it is 
not expected to be well described by a jellium-like model. 

A face-dependent approach extension of the stabilized- 
jellium model consists in obtaining the self-consistent 
electron density by adding to the constant poteptiial 
((5w)^_5 a structure-dependent corrugation factor£3't£l 
This procedure yields an increased surface energy [hori- 
zontal dashed-dotted lines of Figs. 2 and 3], which in the 
case of Al is found to be close to the experimental result. 

Figs. 4 and 5 exhibit our calculated stabilized-jellium 
work functions for slabs of Al and Li, respectively, as a 
function of the thickness L of the slab, together with first- 
principles thin- film calculations. As in the case of the sur- 
face energy, a procedure similar to that of Eq. ( |l5|) yields 
a work function [represented by horizontal solid lines] 
that agrees within less than 0.1% with the result we also 
obtain after solving Eq. (^ for the semi-infinite medium, 
a precision that is difficult to achieve by a fitting proce- 
dure. For L ~ 0.5Af, the QSE yields oscillations with 
relative amplitudes of ^ 20% and ~ 10% for Al and Li, 
respectively. For Al both the amplitude and the oscilla- 
tion pattern are comparable to those exhibited by atom- 
istic calculations. In the case of a 3-layer film of Al(lll), 
the slab width is L 4(Af/2). Hence, the stabilized- 
jellium model predicts a minimum for this film, which is 
in reasonable agreement with the deep minimum exhib- 
ited by atomistic calculations with v = Z. In the case 
of Li(OOOl), the stabilized-jellium model predicts a min- 
imum for a one-layer film \L '-^ \ (Af/2)], also in agree- 
ment with the minimum exhibited by first-principles cal- 
culations with V = \. Finally, we note that adding a 
structure-dependent corrugation factor to the slabilized- 
jellium (<5w)j^g constant potential yields a smaller value 
of the work function [horizontal dashed-dotted lines of 
Figs. 4 and 5], which in the case of Al is in reason- 
able agreement with the experiment. For Li, both the 
stabilized-jellium model and first-principles calculations 
predict work functions that are well above the experi- 
mental result. 

For given values of the equilibrium-density parameter 
Ts and the valence Z, all these calculations have been 
carried out with the core radius [characteristic of each 



metal] that is obtained from the bulk stability condition 
expressed by Eq. M). However, while at the equilibrium 
density n of Eq. (S) the infinite homogeneous system is 
stable, at this density a finite system is not stable against 
changes of the background density, i.e.. 



d{E/N) 
dr. 



7^0, 



(16) 



where N represents the particle number. Instead, there is 
a modified equilibrium-density parameter r*, which sta- 
bilizes the finite system. This modified parameter de- 
pends on the size L of our system and is expected to 
approach as L — > oo. 

Fig. 6 shows the result of our full self-consistent Kohn- 
Sham calculations of the deviation function 
of the thickness L of the slab. These calculations indi- 
cate that there is a self-compression effect, which is more 
pronounced when the two surfaces are separated by a 
multiple of ^ Af/2. 

The self-compression effect exhibited in Fig. 6 may be 
approximately predicted with use of the LDM, a sis 
model to evaluate the total energy of a finite systemES 
In this model, the energy is the sum of a volume term 
(the bulk energy per unit volume, ncbuik, times the vol- 
ume) and a surface term (the surface energy, a, times the 
transversal area): 



E 



LDM 



nebulk V + (tA. 



(17) 



For fixed rc, and evaluated at the bulk equilibrium- 
density parameter r^, 



djELDAl/N) 

dr. 



A da 



-.^>0. 
dr. 



(18) 



The first term is positive, as can be found from the data 
in Table I of Ref. For a fixed slab width i, the 
second term is also positive, and the surface term self- 
compresses, therefore, stabilized-jellium slabs. The devi- 
ation of the electron density parameter r* obtained from 
the LDM stability condition 



d{ELDM/N) 

dr. 



= 



(19) 



with respect to the hulk equilibrium density parameter 
rs is also plotted in Fig. 6, showing that the LDM pro- 
vides a nice average of our self-consistent Kohn-Sham 
calculatiems, as previously demonstrated in the case of 
clusters. c2l 

In Ref. ^, thin films of Be with 1-3 layers were exam- 
ined and a jellium version of a crystalline calculation was 
considered. The electron density parameter r* needed to 
define each slab was derived from the optimized (relaxed) 
structural parameters. The results reported in Ref. ^ are 
in agreement with the compression effect we report here, 
with r* increasing with the number of layers and ap- 
proaching the hulk equilibrium-density parameter as 
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L — ^ cxo. These results show deviations of the electron 
density parameter, r* — Vg of ^ 3.2%, 1.9% and 0.9% 
for thin films with 1, 2 and 3 layers, respectively. This 
is in agreement with our stabilized-jellium calculations, 
which in the case of thin films with ^ 2 layers of Li and 
Al predict (see Fig. 6) differences between r* and of 
~ 1.6% and 1.8%, respectively. The self-compression of 
structural parameters in ultra-thin crystalline films has 
also been discussed in terms of the so-called coordination 
model which, however, seems to fail in some cases (see, 
e.g., Ref. |lD. 

Finally, we note that if for each value of L the corre- 
sponding equilibrium-density parameter r* is taken, in- 
stead of the bulk parameter , modified surface energies 
and work functions are obtained which are quite similar 
to those displayed in Figs. 2-5. This is in contrast with 
the discussion of Ref. 0. 



IV. CONCLUSIONS 

We have modeled thin films of two simple metals, alu- 
minum and lithium, using the stabilized-jellium model, 
and have studied the convergence of some physical quan- 
tities (work function and surface energy) to the semi- 
infinite planar-surface results. We have found the same 
oscillatory behavior which is typical of the QSE in jel- 
lium. Although this behavior also shows up in atomistic 
first-principles thin-film calculations, the clean QSE of 
continuous background models is obscured in the more 
realistic calculations. A trend consisting in surface en- 
ergy minima coinciding with work function maximaj-w&Si 
reported for first-principles crystalline calculations J13a 
However, within the stabilized-jellium model we have 
found minima and maxima of both quantities at the same 
positions [as also reported in Ref. || from first-principles 
for Al(lll)]. On the other hand, we have found that both 
the absolute and the relative amplitude of stabilized- 
jellium QSE oscillations are larger for aluminum than for 
lithium, in agreement with first-principles evaluations. 
The disagreement between our stabilized-jellium results 
for lithium and the more realistic atomistic all- electron 
calculations cannot be attributed to some property of the 
pseudopotential, and simply shows that this metal does 
not display free-electron behaviour. 

Stabilized-jellium slabs of aluminum and lithium have 
been found not to be stable at the bulk equilibrium den- 
sity, the size-dependent equilibrium density being larger. 
This self-compression effect, which was already known to 
exist for clusters, has been found to become more im- 
portant as the slab width decreases. Both LDM and 
full self-consistent DFT calculations have shown a larger 
self-compression for aluminum than for lithium, which 
is a consequence of the larger surface energy of the for- 
mer material. The self-compression of thin simple-metal 
films is a general rule that is also exhibited by atomistic 
first-principles calculations, where the unitary cell of thin 



films is found to be slightly smaller than that of the bulk 
solid. 

The stabilized jellium model is computationally as sim- 
ple as the jellium model; however, for the two high- 
density metals we have considered, it is much more real- 
istic. In particular, we have found it to be more realistic 
for aluminum than for lithium. The stabilized-jellium 
model is adequate to obtain general qualitative conclu- 
sions and an understanding of trends of simple metals 
but, obviously, is unable to provide precise quantitative 
conclusions on particular metals. These can only be ex- 
tracted from the now standard first-principles, but com- 
putationally more demanding, calculations. 
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FIG. 1. Normalized valence-electron density in the jellium 
model (solid line) and in the stabilized-jellium model (dashed 
line) for a slab of Al (r^ = 2.07) with thickness L = 2\f- The 
background density is represented by the dark area. The fig- 
ure also displays the effective potential 14/ / {z) in each model 
(solid line for the jellium model and dashed line for the sta- 
bilized-jellium model). 



FIG. 2. Surface energy and QSE in aluminum (rs = 2.07). 
Large vertical marks across the horizontal axis show the 
widths of unrelaxed fee Al(lll) slabs with 1^ = 1, 12 atomic 
planes. The width L is given by L = !/(\/3/3)a, a being 
the lattice parameter a = (16 tt Z/3)^^'^rs. The solid oscil- 
lating line shows our calculated surface energy of fiat stabi- 
lized-jellium slabs. Solid and dashed-dotted horizontal lines 
represent our calculated surface energy of semi-infinite fiat Al 
(solid line) and fee Al(lll) (dashed-dotted line) stabilized jel- 
lia. The zero-temperature extrapolation of the experirueatal 
liquid-metal surface tension of Ref. ^ divided by I.2JHJ is 
represented by an horizontal arrow. For comparison, atom- 
istic first-principles calculations from Refs. ^ and|^ are also 
displayed, by solid circles and triangles, respectively. The sur- 
face energies of Ref. ^ were obtained using the self-consistent 
pseudopotential method combined with an independent cal- 
culation of the bulk energy per electron. The surface ener- 
gies of Ref. ^ were obtained within an all-electron scheme 
with the use of a linear-combination-of-gaussian-type-orbitals 
fitting function (LCGTO-FF) and with the bulk energy per 
electron extracted from the slab calculations. Dashed lines 
are to guide the eye. 



FIG. 3. Surface energy and QSE in lithium (r^ = 3.24). 
Large vertical marks across the horizontal axis show the 
widths of unrelaxed hcpr^ti(OOOl) slabs with v = 1, 12 
atomic planes [c/a = 1.6463, which corresponds to Ts — 3.13]. 
The slab width is L = va/2 and the structural-parameters' 
ratio c/a = [iQ Z {vs/a)^. The solid oscillating line 
shows our calculated surface energy of flat stabilized-jellium 
slabs. Solid and dashed-dotted horizontal lines represent our 
calculated surface energy of semi-infinite fiat Li (solid line) 
an hep Li(OOOl) (dashed-dotted line) stabilized jellia. The 
horizontal arrow has the same meaning as in Fig. 2. For 
comparison, atomistic all-electron calculations from Refs. ^ 
and B9 are also displayed, by solid circles. These surface ener- 
gies were obtained with the use of a LCGTO-FF and with the 
bulk energy per electron extracted from the slab calculations. 
Dashed lines are to guide the eye. 



FIG. 4. Work function and QSE in aluminum (rs = 2.07). 
All symbols have the same meaning as in Fig. 2. For com- 
parison, atomistic all-electron calculations from Refs. ^ and 
[451 are also displayed, by solid squares and rhombs, respec- 
tively. The work functions of Ref. Biwere obtained within the 
LCAO scheme, and those of Ref. 43 were obtained with the 
use of surface linearized augmented plane waves (SLAPW). 
The experimental polycrystalline work function of Ref. ^ is 
represented by an horizontal arrow. 
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FIG. 5. Work function and QSE in lithium {r, = 3.24). All 
symbols have the same meaning as in Fig. 3. The experimen- 
tal polycrystalline work function of Ref. ^ is represented by 
an horizontal arrow. 

FIG. 6. Relative difference between the actual equilib- 
rium-density parameter r* and the bulk density parameter 
Vs for aluminum (dashed lines) and lithium (solid lines) sta- 
bihzed-jellium films, as a function of the slab width L. 



